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Abstract 

The  3-stage  magnetic  pulse  compressor  MAG  1-D  has  been  mod¬ 
eled  with  a  code  which  calculates  the  hysteresis  behavior  of  magnetic 
cores  placed  within  a  circuit.  The  voltage  waveforms  are  shown  to 
compare  favorably  with  data  taken  by  Barrett.  A  speculative  modi¬ 
fication  to  MAG  1-D,  suggested  by  the  reset  behavior  of  the  model, 
is  discussed. 


Introduction 

This  paper  is  devoted  to  numerical  simulations  of  the  magnetic 
pulse  compressor  MAG  1-D,  which  is  used  in  the  linear  electron  ac¬ 
celerator  ETA-II.  A  circuit  diagram  of  MAG  1-D  is  shown  in  Fig.  1, 
and  the  parameter  values  are  listed  in  Table  1.  Very  schematically, 
the  device  operates  as  follows.  The  energy  storage  capacitor  Co  is 
initially  charged  to  25.5  kV.  After  the  thyratron  fires,  the  energy  is 
transferred  to  the  precompression  capacitor  Cj,  with  a  charge  time 
of  about  3.8  /ts.  The  magnetic  switch  L\  then  saturates,  allowing 
the  energy  to  flow  through  the  1:10  transformer  and  to  charge  the 
capacitor  Ci  (charge  time  about  1  /ts).  After  the  switch  £2  satu¬ 
rates,  the  energy  is  transferred  to  a  coaxial  cable  pulse-forming  line 
(electrical  length  35  ns),  which  has  a  charge  time  of  about  200  ns. 
At  this  point,  the  third  magnetic  switch  saturates  and  the  energy  is 
applied  to  the  load,  which  is  here  treated  as  a  simple  resistor.  The 
unabsorbed  energy  then  reflects  back  through  the  circuit.  The  cores 
in  MAG  1-D  are  reset  by  two  circuits,  as  shown  in  Fig.  1.  All  cores 
are  wound  with  Metglas  2605CO. 

We  model  this  device  making  use  of  a  core  model  discussed  ear¬ 
lier  (1).  Magnetic  hysteresis  is  treated  by  a  rate-dependent  theory 
developed  recently  by  Hodgdon  [2].  Rate  dependence  is  found  to 
play  an  important  role  in  the  second  stage  and  a  crucial  role  in  the 
third  stage. 

The  voltage  waveforms  are  compared  with  data  taken  by  Bar¬ 
rett  [3].  At  the  beginning,  there  is  essentially  complete  agreement, 
while  at  the  end,  the  main  pulse  is  accurate  to  within  about  20%. 
Only  a  single  cycle  of  duration  200  /is  was  simulated. 

The  first  reset  current  is  found  to  undergo  wide  swings  with  a 
period  of  about  70  /is  and  an  amplitude  of  about  100  amps.  Con¬ 
sequently,  the  first  stage  switch  does  not  reset  properly,  and  the 
precompression  capacitor  (input  to  the  first  stage)  undergoes  multi¬ 
ple  voltage  oscillations.  To  investigate  whether  this  reset  behavior 
could  be  corrected,  the  model  was  modified  so  as  to  contain  a  single 
isolated  reset  rod.  The  reset  current  then  settled  down,  and  the  pre¬ 
compression  capacitor  showed  only  a  single  late-time  oscillation.  The 
behavior  of  the  pulse  itself  (for  one  cycle)  was  practically  unaffected. 


Model 


In  the  core  model  [1]  employed  here,  the  magnetic  field  within 
the  cores,  in  both  the  switches  and  the  transformer,  is  taken  to  be 
azimuthal  and  to  depend  only  on  time,  and  the  displacement  current 
is  assumed  negligible.  Thus  the  voltage  across  the  core  can  be  shown 
to  satisfy  U(t)  =  £(t)  d,I{t)/dt,  where  the  time-dependent  inductance 
is 


£(f)  = 


Kh  n(r,t)  dr  w 

2*  Jr,„  r 


N  h  .r out .  .  , 
— in(  — )/i(t). 
I*  rin 


(1) 


Here  /» is  the  time-dependent  permeability,  h  is  the  height  of  a  core, 
Nw  is  the  number  of  windings,  and  r<n  and  rout  are  the  inner  radius 
and  outer  radius,  respectively.  We  have  ignored  the  radial  depen¬ 
dence  of  the  fields,  which  typically  is  an  effect  of  order  10%.  The 
magnetic  field  itself  evolves  according  to  the  derivative  of  Ampere’s 
law,  or  B  =  nNV/(2xfL),  with  fi  =  (rou(  -  rin)/ hi(rmt/rin). 

The  complications  of  hysteresis  reside  in  /t,  which,  in  the  Hodgdon 
theory,  depends  on  B,  B,  and  H  according  to  an  equation  of  the  form 


ti(B,H,B)  =  {asign(B)  [/(B)  -H]  +  g(B,B)}\  (2) 

where  a  is  a  constant  and  /  and  g  are  functions  particular  to  the 
magnetic  material.  The  fact  that  the  permeability  depends  on  the 
sign  of  B  leads  to  hysteresis.  This  theory  reproduces  many  of  the 
known  features  of  hysteresis  phenonmena,  such  as  loop  closure  near 
saturation,  the  formation  of  minor  loops,  and  loop  broadening  with 
increasing  B.  The  functions  appropriate  to  Metglas  2605CO  are 
given  in  the  Appendix.  In  a  previous  application  [1],  this  combi¬ 
nation  of  core  model  and  hysteresis  model  was  shown  to  give  quite 
satisfactory  agreement  with  experiment. 

We  employ  a  simple  thyratron  model,  which,  for  positive  currents, 
is  a  negative  500  v  source.  This  accounts  for  the  forward  drop  across 
the  thyratron.  For  negative  currents,  the  thyratron  is  modeled  as  a 
1000  ohm  resistor  in  series  with  the  same  source. 

The  coaxial  cable  PFL  is  treated  as  a  series  of  N  stages,  each 
containing  a  parallel  capacitor  C/N  and  resistor  NG,  with  a  series 
inductor  L/2N  and  resistor  R/2N  on  each  side  (cf.  Table  1).  Here 
R  accounts  for  the  total  ESR  of  the  line,  and  G  accounts  for  the 
conductivity  of  water.  The  simulations  reported  here  have  N  = 
10  stages. 

The  equations  of  the  model  are  set  up  as  a  system  of  ODEs.  In  the 
problem  at  hand,  the  system  has  54  variables,  including  magnetic 
fields  within  the  cores.  Owing  to  the  presence  of  rate-dependent 
hysteresis,  the  equations  have  extreme  nonlinearities.  The  system  is 
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Results 

We  now  turn  to  a  full  200  /is  simulation  of  the  device.  At  t  =  0, 
the  intermediate  energy  storage  capacitor  Co  is  charged  to  25.5  kV, 
and  the  reset  currents  are  assigned  their  static  values.  Figure  2 
shows  the  calculated  voltage  across  C\,  along  with  data  taken  by 
Barrett.  The  main  pulse  is  very  close  to  the  data,  but  the  reflected 
pulse  is  about  20%  smaller  than  experiment.  In  Fig.  3  this  voltage 
is  shown  on  a  long  time  scale.  The  most  conspicuous  feature  here 
is  the  sudden  reversals  at  35,  100,  and  175  ps,  which  are  caused  by 
resaturations  of  the  magnetic  switch  in  this  stage.  Similar  reversals 
are  seen  experimentally. 

Figure  4  shows  the  B-H  loop  of  the  first  magnetic  switch.  Note 
that  the  curve  balloons  to  the  right,  tucks  in  somewhat,  and  then 
saturates.  After  the  second  saturation,  however,  the  fields  wander 
around  and  fail  to  reset  properly.  The  reason  for  this  can  be  seen 
in  Fig.  5,  which  shows  the  reset  currents  as  functions  of  time.  The 


Fig.  2.  Voltage  across  Ci  (short  time  scale).  The  dotted  line  indi¬ 
cates  Barrett’s  data  [Ref.  3]. 

Table  1.  MAG  1-D  parameter  values 


Stage  1  Stage  2  Stage  3 


Number  of  cores 

3 

4 

3 

Number  of  windings 

2 

3 

1 

Inner  radius  (cm) 

8.89  8.73 

11.43 

Outer  radius  (cm) 

15.24  18.73 

18.73 

Core  height  (cm) 

5.08  5.08 

5.08 

Saturated  inductance  (nH) 

141 

631 

27.3 

Stray  inductance  (nH) 

50 

240 

5 

Ri  (mfl) 

.375  6.19 

.77 

Packing  factor 

.67 

.67 

.55 

Transformer  parameters: 

PFL  parameters: 

Number  of  turns:  10 

Inner  radius:  16.19 

cm 

Number  of  series  cores:  4 

Outer  radius:  22.23  cm 

Core  height:  5.08  cm 

Length:  116.8  cm 

Inner  radius:  8.89  cm 

C  =  17.3  nF 

Outer  radius:  16.35  cm 

L  =  74  nH 

Saturated  inductance:  52.9  nH 

R  =  .07  fl 

Packing  factor:  .67 

G  =  7770  fl 

Primary  ESR:  .3  mfl 

Secondary  ESR:  89  mfl 

Lo  =  1.5  #iH  Rl0  =  4.35  mfl 

V{*  =  20  v 

Co  =  2  /XF  Rc0  =  21  mfl 

R\'  =  in 

Ci  =  2.1  nF  RCl  =  7  mfl 

L[e  =  150/rH 

C2  =  19.3  nF  Rl  =  740  mfl 

Fjre  =  50  v 

-filed  =  2.1  fl  R2C7  =  46  mfl 

R S'  =  2.5  n 

£io.d  =  10  #*H 

L 5'  =  4.35  mH 

first  reset  current  oscillates  greatly  and  reaches  negative  maxima 
in  synchronization  with  the  current  through  the  magnetic  switch. 
This  problem  seems  not  to  be  present  in  the  second  reset  current.  It 
reaches  a  maximum  about  when  the  third  magnetic  switch  saturates, 
and  then  quickly  settles  down. 

Figure  6  shows  the  voltage  across  the  water  capacitor  C2,  along 
with  the  experimental  trace.  It  has  a  charge  time  of  1.2  /is,  a  peak  of 


Time  (microsec) 


Fig.  3.  Voltage  across  C\  (long  time  scale). 


Fig.  4.  B-H  loop  of  magnetic  switch  in  stage  1. 


Fig.  5.  First  reset  current  (full  line)  and  second  reset  current 
(dashed  line). 


230  kV,  and  a  discharge  time  of  about  200  ns,  all  in  reasonable  agree¬ 
ment  with  Barrett’s  measurements.  The  B-H  loop  of  the  magnetic 
switch  in  this  stage  is  shown  in  Fig.  7.  Before  the  first  saturation,  it 
balloons  about  twice  as  far  to  the  right  as  the  previous  switch,  as  it 
experiences  more  severe  rate  dependence.  Afterwards,  it  resaturates 
and  then  appears  to  reset  normally. 

The  voltage  across  the  output  of  the  PFL  is  shown  in  Fig.  8.  Its 
maximum  (210  kV)  is  about  10%  higher  than  experiment.  Note  that 
it  falls  in  about  10  ns  to  half  maximum,  momentarily  reverses,  find 
then  falls  in  another  10  ns  to  zero.  This  behavior  can  probably  be 
traced  to  the  interplay  of  forward  and  backward  waves.  A  similar 
pattern  is  observed  in  experiment.  The  multiple  reflections  following 
the  pulse  have  a  period  of  twice  the  electrical  length  of  the  PFL,  or 
70  ns.  These  also  show  up  on  the  experimental  trace. 

Figure  9  shows  the  B-H  loop  of  the  magnetic  switch  in  the  third 
stage.  It  is  wider  by  a  factor  of  four  than  that  in  the  previous  switch. 
The  multiple  reflections  in  the  PFL  show  up  here  in  a  picturesque 
way,  and  the  switch  appears  to  reset  properly. 

Figure  10  shows  the  calculated  and  observed  voltages  across  the 
load.  The  calculated  peak  of  118  kV  is  about  20%  higher  than  ex¬ 
periment,  but  the  widths  (about  60  ns)  agree  well. 

Turning  to  the  final  energy  balance,  we  find  that  the  load  dissi¬ 
pates  467  joules.  The  greatest  unwanted  energy  loss  occurs  in  the 
thyratron,  which  dissipates  34  joules.  Hysteresis  losses  in  the  third 
stage  account  for  the  next  most  serious  loss  (29  joules). 


Fig.  6.  Voltage  across  C2.  The  dotted  line  indicates  Barrett’s  data. 


Fig.  7.  B-H  loop  of  magnetic  switch  in  stage  2. 


Fig.  8.  Voltage  across  output  of  PFL.  The  dotted  line  indicates 
Barrett’s  data. 


^  Y  7  0  ~  ~  ^ 
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Fig.  9.  B-H  loop  of  magnetic  switch  in  stage  3. 


Fig.  10.  Voltage  across  load.  The  dotted  line  indicates  Barrett’s 
data. 
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Conclusions 


Time  (microsec) 

Fig.  12.  Voltage  across  C\  (alternate  case). 


These  simulations  show  that  the  MAG  1-D  modulator  can  be  mod¬ 
eled  with  quite  reasonable  accuracy.  Hence  the  code  could  serve  as 
a  useful  guide  in  trouble-shooting  and  in  the  design  of  possible  mod¬ 
ifications  to  the  device. 

As  a  specific  example  of  such  interplay  between  model  and  experi¬ 
ment,  the  results  indicate  that  the  first  stage  does  not  reset  properly 
in  the  base  case.  An  alternative  reset  mechanism  was  set  up  and 
shown  to  give  better  performance. 
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Appendix 

In  the  version  of  the  Hodgdon  theory  used  here,  the  material  func¬ 
tions  /  and  g  have  the  form 

(  Ai(tan(A2B))p  -f-  /e*B,  0  <  B  <  Bbp\ 

/(B)  =  <  Ai(tan(A2B(,p))p  +  fexBbp  +  (B  -  Bbp)/p„  B  >  Bbp\ 

{  -/(-B),  B  <  0; 

(Al) 

,(*,*)  =  |  «*>['  -  *  "I*'  <  3 

[  f'(B),  otherwise. 

m 

The  function  c(B),  which  controls  the  rate  dependence,  is  given  by 


fl  +  ci|B|,  |B|<Bi; 

c(B)  =  {  1  +  c,B!  +  c2(|B|  -  Bi),  Bi  <  |B|  <  B2; 

{  1  +  c,Bi  +  c2(B2  -  BO  +  c3(|B|  -  B2),  |B|  >  B2. 

(A3) 

For  Metglas  2605CO,  the  values  of  the  constants  are,  in  mks  units: 

Bc|  =  1.7  A!  =  .05809  cx  =  1.1  x  10~4  Bt  =  1  x  106 

Bbp  =  1.75  A2  =  .8941  c2  =  6.17  x  10"'  B2  =  3  x  106 

/„  =  .7958  A3  =  -26.23  c3  =  3.061  x  10~6  B3  =  2.1  x  107 

p  =  2  A4  =  .3101  p,  —  .582po  a  =  10 


Alternate  Case 

To  see  whether  the  reset  problem  observed  in  the  previous  run 
could  be  improved,  at  least  conceptually,  the  modulator  was  imag¬ 
ined  to  be  reconfigured  such  that  an  isolated  reset  rod  passed  axially 
through  the  centers  of  all  the  cores.  (This  is  possible  only  if  all  stages 
are  configured  with  their  fluxes  pointing  in  the  same  direction.)  In 
the  simulation  the  reset  rod  was  connected  to  a  voltage  source  of 
20  v,  with  a  1  ohm  resistor  and  a  large  inductor  (10  mH)  in  series. 
All  other  run  conditions  remained  the  same  as  above. 

The  resulting  reset  current  is  stable  and  well-behaved,  rising  by 
50%  and  then  settling  back  toward  its  equilibrium  value  of  20  amps. 
Hence  the  first  switch  also  settles  down,  as  shown  by  its  B-H  loop 
in  Fig.  11.  Note  that  the  cores  have  not  quite  fully  resaturated  even 
after  200  ps;  some  further  refinement  of  the  reset  parameters  might 
be  necessary  to  achieve  a  complete  reset.  The  voltage  across  Ci  now 
undergoes  only  a  single  late-time  oscillation,  at  about  100  ps,  as  seen 
in  Fig.  12.  This  figure  should  be  compared  with  Fig.  3,  showing  the 
multiple  oscillations  observed  with  the  original  reset  scheme.  The 
other  switches  and  the  PFL  perform  much  as  in  the  previous  case. 
The  voltage  across  the  load  has  virtually  the  same  shape  as  in  the 
previous  case. 
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